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Nuclear magnetic resonance (NMR) has been widely used as a 
demonstrative medium for showcasing the abiUty for quantum 
computations to outperform classical ones. A large number of 
such experiments performed have been implementations of the 
Deutsch- Jozsa algorithm. It is known, however, that in some cases 
the Deutsch-Jozsa problem can be solved classically using as 
many queries to the black-box as in the quantum solution. In this 
paper we describe experiments in which we take the contrasting 
approach of using NMR as a classical computing medium, treating 
the nuclear spin vectors classically and utilising an alternative 
embedding of bits into the physical medium. This allows us to 
determine the actual Boolean function computed by the black-box 
for the n = 1,2 cases, as opposed to only the nature (balanced or 
constant) as conventional quantum algorithms do. Discussion of 
these experiments leads to some clarification of the complications 
surrounding the comparison of different quantum algorithms, 
particularly black-box type algorithms. 
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1 INTRODUCTION 



Nuclear magnetic resonance (NMR) experiments are generally conducted 
using bulk samples and hence the manipulating radiofrequency (rf) pulses 
and the detection signal have to be regarded in the context of an ensemble 
average of the underlying nuclear magnetic spin dynamics. Theoretically, 
this situation is successfully dealt with by a density matrix approach. Since, 
however, the idea of quantum computation is based on the concept of being 
able to manipulate the spin dynamics on the basis of pure quantum spin 
states, there have been various attempts at implementing quantum computation 
algorithms using the experimental conditions and restrictions of NMR by 
adopting pseudo-pure spin state based approaches. 

The most commonly implemented quantum algorithm, both in NMR and 
in general, is the one due to Deutsch and Josza [S, 23, 9, 10, 14, 20, 26, 3, 24, 
19, 21, 11, 29, 25, 16, 15, 18]. The various NMR implementations differ by: 
the underlying spin quantum numbers 5(S' = l/2or5> 1/2); the initial 
spin states (thermal equilibrium state or pseudo pure state); the algorithmic 
implementation of the problem (Collins [ ] or Cleve [ ] defining the number 
of qubits necessary to operate a given DJ problem size). 

In the context of computation, NMR has in the past been exclusively used 
for implementing quantum computations. However, it also has potential as a 
classical computing medium, where the rich state space can be fully utilised to 
perform classical operations [28]. The Deutsch- Jozsa problem, long touted as 
a simple and key example of the ability of quantum computing to outperform 
classical computing, has more recently been shown to be de-quantisable in 
some cases — i.e., efficient classical solutions can be formulated [1, 6]. 

In this paper, we describe the implementation of the n — 1 and n = 2 
de-quantised solutions for the DJ problem in a classical NMR computation. 
The process of implementing this solution highlights key aspects of quantum 
algorithms and computation, and we discuss these in detail. In particular, we 
emphasise the separation between three nested 'layers' of any quantum algo- 
rithm: the problem formulation, the algorithm formulation, and the physical 
implementation. In general these levels are independent, but certain conditions 
on the relationship between levels must be satisfied. Specifically, a partic- 
ular algorithm applies only to a specific problem formulation, and for each 
algorithm a choice of embedding into the physical medium must be made 
in order to implement it. Further, for 'oracle' or 'black-box' problems such 
as the Deutsch-Jozsa problem the comparison of different formulations of 
the problem requires discussion of the ability to embed the black-box from 
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one formulation into the other. All these issues are subtle and require further 
discussion. 

2 PROBLEM FORMULATION 

The standard formulation of the Deutsch-Jozsa problem [' ^] is as follows. 
Let / : {0, 1}" — > {0, 1}, and suppose we are given a black-box computing 
/ with the guarantee that / is either constant (i.e. for all xi,X2 E {0, 1}" 
we have f{xi) — f{x2)) or balanced (i.e. f{x) = for exactly half of all 
possible inputs x £ {0, 1}"). The problem is to determine, in as few black-box 
calls as possible, whether / is constant or balanced. The obvious classical 
algorithm must examine one more than half the input bit-strings and thus 
requires 2"^^ + 1 black-box calls, while the quantum solution requires only 
one. 

There is, however, an important subtle issue: the classical and quantum 
problems are slightly different. In one case we are given a classical black- 
box Cf computing /, and in the other we are given a unitary black-box Uf, 
operating in a 2"+^ dimensional Hilbert space H2"+i > computing Uf \x) \y) = 
\x) \y ® f{x)) where '0' denotes addition modulo 2. As such, it is bending the 
truth a little to say that 'the problem' can be solved more efficiently quantum 
mechanically than classically. 

The possibility in the n = 1, 2 cases to de-quantise the quantum solution 
to give an equally good classical algorithm [ , ] results from working with 
a formulation of the problem in which we are given a higher dimensional 
classical black-box. In this case, complex numbers are used as 'complex 
bits' — a classical analogue of a qubit — and the black-box Cf operates (in the 
n — 2 case) as follows: 

/ai _ [oi + (-l)/(oo)e/(io)5^,] x 

U + b2i)~\ a2 + (-l)/(io)®/(")62* J ■ 

(1) 

In general, a particular algorithm, be it classical or quantum, solves a particular 
formulation of the problem; i.e., it determines if a black-box of a specific type 
which computes in some reasonable form f is balanced or constant. Hence, 
it seems at least some of the apparent difference in powers of the classical 
and quantum solutions comes from the slightly different formulation of the 
problems, i.e., the different 'powers' of the black-boxes. 

A comment should be made about what it means to compute / in some 
reasonable form. Since quantum computing requires unitarity, the simplest 
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and perhaps most natural way to compute / is with an /-controlled-NOT gate; 
indeed this was the original method proposed by Deutsch [ ]. One thus needs 
to be careful of claims that the separate output qubit for Uf is not needed [8] 
and that the alternative quantum black-box U'^ \x) = (—1)-'''^^) |a;) can equally 
well be used. Rather, it seems that U'j does not reasonably compute / as 
absolute phase factors have no physical significance and it is hence impossible 
to characterise which boolean function / the black-box 'computes' by trying 
various inputs. The same issue is not present, however, in the classical de- 
quantised solution because phase factors are measurable in this case, and thus 
C f can be seen to compute /, albeit in a slightly non-standard way. 

2.1 Black-box embeddings 

The quantum black-box (represented by Uf) is often considered an embedding 
of the classical black-box computing / [ ]; if this were true it would be more 
reasonable to view the quantum solution as solving the original problem. This, 
however, is a misunderstanding which helps hide the differences between the 
classical and quantum formulations of the problem. For this to be an embed- 
ding we would require a map e : {0, 1}" ^ ^g(2i+i) (where denotes 
an embedding) which preserves the structure of the computed function. That 
is, for X e {0, 1}", e{f{x)) = Uf{e{x)). In other words, the computational 
states we assign to represent the binary bits and 1 must behave as expected 
under Uf given that Uf should compute the function /. 

However, the requirement of the unitarity of Uf makes such an embedding 
impossible. This can be verified by considering any constant boolean /: such 
a function is not bijective, so no bijective Uf preserving the required structure 
can exist. Because no unitary embedding is possible, it seems more suitable 
to consider the quantum solution as a method to solve an analogue of the 
classical problem, rather than a more efficient solution to the classical problem. 
One is forced to conclude that the typical claims comparing the quantum and 
classical solution are, in fact, not valid. Any comparison of the problems 
should take into account the differences in complexity of the black-boxes [2]. 

Interestingly, there is an embedding between the quantum black-box Uf 
and the de-quantised black-box C/, so it is not as unreasonable to compare 
the solutions using these black-boxes as it is to compare the quantum and one- 
dimensional classical solutions. By realising that the original classical solution 
and the de-quantised classical solution are not solving the same problem, we 
see there is no explicit contradiction with claims that 2"^^ + 1 black-box calls 
is the best that can be done in the original classical problem [27]. 
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3 ALGORITHM FORMULATION 



Given a particular formulation of the problem, the algorithm formulation 
involves determining the input for the black-box, and what operations are 
required to determine the nature of / from the output of the black-box. 

In the standard quantum solution [7, 13], we initially prepare our system 
in the state |00) and then operate on it with a three-qubit Hadamard gate, 
H^^, to get: 

1 



i7«^|00) |1) = 



2 S 



X - = 



(2) 



After applying the /-controlled-NOT gate Uf we have 



Uf 



-2 E 



N>i->= E 

x£{()A}^ 



x)\-). 



(3) 



By applying a final 3-qubit Hadamard gate to project this state onto the 
computational basis we obtain the state 



(_i)/(oo) ij(oo) e /(lo)) ® 1/(10) ® /(II)) |i) 



(4) 



If we measure both the first and second qubits we can determine the nature of 
/: if both qubits are measured as 0, then / is constant, otherwise / is balanced. 
This result is correct with probability one. 

The de-quantised solution works in a similar way, but uses complex num- 
bers as two-dimensional complex bits. Using the black-box Cf defined previ- 
ously, the algorithm proceeds as follows. We set z — zi — Z2 — I + i, apply 
Cf and multiply by z to project onto the computational basis to obtain the 
result: 



'(-i)"™='\ A-i)"»»'A ,,,, , , 

= if / IS constant. 



'(_1)/(00)^A /(_l)/(00) 



X < 



zz , 

'(_l)/(oo)^2^ 

zz I 




if / is balanced. 



(5) 
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FIGURE 1 

Creating and observing a NIVIR signal. After reacliing thermal equilibrium magnetisa- 
tion in the external magnetic field, the magnetisation vector S is flipped to the sy-plane 
of the laboratory frame of reference by applying a radio-frequency (rf) pulse of suitable 
duration. After the rf pulse, the NMR signal is detected in the xy-plane in the form of 
a time-domain signal, the free induction decay (FID) which is recorded by a receiver 
coil, here assumed to be placed in the a;-direction. The FID is converted into the 
frequency-domain spectrum by a Fourier transform. 



By checking both of the resulting complex numbers, we can determine 
whether / is balanced or constant with certainty. If both complex numbers 
are imaginary then / is constant, otherwise it is balanced. In fact, the ability 
to determine if the output numbers are negative or positive allows us to deter- 
mine the value of /(OO) and thus which Boolean function / is; the quantum 
algorithm is incapable of doing this [27]. 

The ability to de-quantise the n = 1 and n — 2 solutions suggests that, at 
least in these cases, the power of the quantum algorithm comes from exploiting 
the two-dimensionality of the black-box, rather than from quantum mechanical 
effects [1]. 

4 IMPLEMENTATION OF CLASSICAL NMR COMPUTING 

For the implementation of the de-quantised algorithm we use real, two- 
dimensional vectors to represent our complex bits (only the direction of the 
vector is of particular importance); this is equivalent to using complex num- 
bers, but more convenient for use in this implementation. In this representation 
our basis bits corresponding to the classical and 1 become (1,0) and (0, 1) 
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respectively. For the implementation we must embed our bits into the physical 
medium; to do so, we use the nuclear magnetisation spin vectors for this 
embedding. Specifically, we embed the complex bits into the xy-plane of the 
rotating frame of reference in a NMR experiment as is illustrated in Fig. I . We 
only consider uncoupled spin species with spin quantum number S = 1/2, and 
we have the advantage that every direction of the magnetisation vector in the 
xy-plane is distinguishable. Working solely with uncoupled spin species, their 
dynamics are fully described by a classical model of magnetisation vectors 
subjected to a range of differnt rotations (pulses) ["''", ']. Hence, at no point 
do we make explicit use of the quantum-mechanical properties of nuclear 
spin systems. In particular, we take (row-vectors represent complex bits, col- 
umn vectors are the nuclear spin vectors; we omit normalisation factors for 
convenience): 

Combinations of the complex bits are naturally taken by the vector addition 
(i.e. the embedding is linear) of the corresponding spin vectors (e.g. (1, 1) — >■ 

= (J))- This mapping is shown in Fig. 2. We emphasise that while the 
same nuclear spin state vectors are used as in conventional quantum computing 
experiments, we choose a different embedding of our algorithm into this spin- 
state vector space which is more flexible for use with complex bits. 

Our sample consists of 99.8% deuterated chloroform with a small amount 
of H2O added. The spins in the residual CHCI3 and H2O in this mixture 
are used for the implementations. The top row in Fig. 3 depicts a conventional 

NMR spectrum of the sample. The implementation for n = 1 only requires 
one spin species to be present. This is most easily achieved by using selective 
excitation pulses, centred at the resonance frequency of the CHCI3 NMR 
resonance. The NMR spectrum obtained by selective excitation, together 
with the corresponding excitation profile of the selective pulses used are shown 
in the middle and bottom traces of Fig. 3. 

4.1 The n — 1 implementation 

The n = 1 implementation relies on the form of the black-box Cf operating as 

Q((ai,6i)) = ((-iy(°)ai,(-l/WM. 

The NMR pulse sequence for implementing this for our sample is shown in 
Fig. 4. The sequence starts with a (selective) tt/2 pulse, flipping the (CHCI3) 
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FIGURE 2 

The embedding from classical two-dimensional bits onto the nuclear magnetisation 
vectors in the xy-plane. Assuming signal detection along the a;-direction, the phases 
of the corresponding NMR spectra are shown, together with the basis bits (0, 1) and 
(1,0) (see text). 



H magnetisation vector from the initial equilibrium z-direction into the xy- 
plane, followed by a sequence of two tt pulses applied to the xy-magnetisation 
vector as required to implement the black-box before the resulting signal is 
detected in the form of a FID. 

The four possible boolean functions for n = 1 may be written /ab where 
A — /(O) and B = /(I). The four corresponding black-boxes can be 
uniformly implemented as the following set of tt pulses, applied to oiu" basis 
bits (see Fig. 2): 



We see that the physical embedding of the black-box fulfils the requirement 
that the function / is reasonably computed due to the ability to distinguish all 
directions of magnetisation vectors in the xy-plane, and thus all valid complex 
bits. In particular, applying the black-boxes for inputs (1, 0) or (0, 1) yields 
((_l)/(o) 0) and (0, (-l)^(^)) respectively as desired. The effect of the four 
black-boxes on the 'basis bits' is shown in Table 1, both in terms of complex 
bits and spin vectors. 

For the actual algorithm, the de-quantised solution can be simplified by 
dropping the use of any equivalent of the final Hadamard operation; there is 
no need to project the result onto the basis states since non-basis states are 
equivalently detectable. In fact, the choice of /-1-45 states as our 'basis bits' 
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FIGURE 3 

NMR spectrum (loo /2tv = -600.13 MHz) of 99.8% deuterated chloroform with a 
small amount of H2O added (top trace). The middle trace shows the NMR spectrum 
after application of a selective 90° pulse, centered around the resonance frequency of 
the CHCI3 NMR signal. The bottom trace displays the excitation profile of the 
selective pulses used. 



(1)90 (41' FID 



FIGURE 4 

The NMR pulse sequence for implementation of the n — 1 case; the flip angles 
effected by a rf pulse are given in fractions of n, the symbols A, B refer to the control 
parameters of the black-box, curved shapes of rf pulses indicate selective pulses [17]. 
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A = /(O) 


B = f{l) 


C/((1,0)) 


C/((0,1)) 
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TABLE 1 

The effect of the four black-boxes on the basis-bit inputs for n = 1. 
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TABLE 2 

The algorithm run with each of the four black-boxes for n = 1. 



(rather than and ly) is convenient because this leaves the system in or ly 
states after the black-box, and these states are particularly easy to distinguish 
(see Fig. 2). 

The effect of the algorithm with each black-box is shown in Table 2. The 
coiTesponding spectra relating to the computation are shown in Fig. 5, along 
with the action of the black-box on the basis-state inputs. In particular for a 
constant / we have I^^ spectra, and for balanced we have I±y spectra. Further, 
because we can distinguish J+x spectra from I^^ spectra (and I^y spectra 
from I-y spectra) we can easily determine which function / the black-box 
represents. 

4.2 The n — 2 implementation 

The idea is that the n — 2 case should be implemented as a natural extension 
of the n = 1 case by expanding it to include two different spin species. The 
function to be implemented is: 

Cji{aMia2,h)) = (ai(-l)/(°"),6i(-l)/(i°)) [a,,b,i-iy^''^^f<^''^) . 
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FIGURE 5 

'^H NMR spectra, CHCI3 resonance, implementing ttie n = 1 case. Tlie top two rows 
show the action of the black-box on the basis-state input. The bottom row shows the 
computation itself, i.e. with (1,1) input. 
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FIGURE 6 

The NMR pulse sequence for implementing the n — 2 case with two different spin 
species, 1 and 2; r'*^^ denotes a pre-acquisition delay of suitable duration. 



With the same mapping of bits as for the n = 1 case, the natural extension of 
the embedding is to have the black-box, now defined by the three parameters 
A = /(OO), B = /(lO) and C = /(lO) /(II), implemented as: 

('^)M5('^)f-45('^)2;-45- 

With this we can see that the first spin is treated as for the n = 1 case, and 
the second spin requires only a single pulse (or none, depending on /). This 
form of the black-box, however, is only valid if one assumes ideal pulses. In 
reality, we now deal with two different spin species with different resonance 
frequencies, hence the receiver can only be on resonance with one species at a 
time. In addition, all rf pulses have non-vanishing durations (in particular the 
rather 'soft' selective pulses with durations of the order of ms) during which 
evolution of magnetisation occurs. Hence, the pulse sequence of the physical 
implementation as shown in Fig. 6 has to account for non-ideal behaviour in 
order to produce an equivalent result to the one obtained from the idealised 
mathematical description. 

The NMR experiment starts with selective pulses centred on spin species 
1 (exactly as for the n = 1 case before), followed by a selective tt/2 pulse 
centred on spin species 2, and a delay of duration t^'^\ If the receiver is 
kept on resonance with species 1 throughout, it is off resonance for species 2. 
Accordingly, after the selective (tx j'l)^ pulse the magnetisation of this species 
will precess in the xy-plane, and by this acquire a phase difference relative to 
the on-resonance signal of species 1 . This precession frequency depends on 
the difference in resonance frequency between the two spin species and hence 
one can easily calculate the correct duration of the delay r'^'-^' that corresponds 
to a dephasing by 37r/2 from the starting 1,^,2 condition to I-y.2- At this point 
in time, the situation is exactly equivalent to the ideal-pulse scenario and the 
data acquisition is started. 

The pulse sequence depicted in Fig. 6 is the simplest form in which the 
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(-1,0)(0, 1) 
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(-1,0)(1,0) 


(-1,0)(0,-1) 
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-i)(i,o) 


(0, 


~1)(0,-1) 



TABLE 3 

The effect of the eight black-boxes on the basis-bit inputs for n = 2. 



required idealised black-box can be carried out in a real NMR experiment. 
For example, one might have started the pulse sequence with a non-selective 
('hard') 7r/2 pulse covering both spin species simultaneously. That would be a 
common preparation step for a spin system, creating a common initial state 
from which to work. However, this would only complicate matters as now 
one would have to take into account dephasing of species 2 magnetisation 
during all selective pulses applied to species 1, and one would have to apply 
the corresponding phase corrections to any pulses applied to species 2 later 
in the sequence. We do not have to start from one common initial state; 
instead our physical implementation avoids all unnecessary such corrections 
and calculations by simply using selective pulses only and preparing species 2 
'just in time' such that only the duration of the pre-acquisition delay needs to 
be calculated. 

As for the n = 1 case, it is easily verified that once again the requirements 
for the black-box to reasonably compute / are satisfied: on a given 'basis' 
input an even number of sign changes indicates a function value of on this 
input, and odd number a value of 1. The results of the black-box on basis 
inputs is shown in Table 3, and the result of the algorithm is shown in Table 4. 
The corresponding input/output spectra are shown, as for the n = 1 case, in 
Fig. 7. 
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FIGURE 7 

NMR spectra, CHCI3 and H2O resonances, implementing the n = 2 case. The top 
four rows show the action of the black-box on the basis-state input. The bottom row 
shows the computation itself, i.e. with (1,1) (1,1) input. 
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TABLE 4 

The algorithm run with each of the eight black-boxes for n = 2. While the idealised 
version of the black-box pulses for species 2 is quoted here, the physically implemented 
pulse sequence gives identical final results. 



5 SUMMARY 

The experiments described in this paper not only show the successful imple- 
mentation of the de-quantised Deutsch-Jozsa algorithm for n — 1,2, but also 
bring attention to some important points regarding quantum computations. 
The experiments confirm that the Deutsch-Jozsa problem can indeed be solved 
classically in some cases using as many black-box calls as the standard quan- 
tum algorithm. During the process of developing our implementation some 
issues about the comparisons of algorithms for different formulations of the 
problem were raised and discussed in detail. We note that, while we explicitly 
demonstrate the ability to solve the n = 1, 2 cases for what is essentially an 
alternative (but equivalent to the quantum) problem formulation, the ability to 
do so for higher values of n reduces to the ability to de-quantise the quantum 
solution for these cases, which is believed not to be possible [ ' ]. 

These experiments reiterate the utility of NMR as a classical computing 
substrate [28]. Specifically, the use of an alternative embedding from the 
problem space onto the uncoupled nuclear spin vectors allowed us to perform 
the classical algorithm and even determine the specific function / computed. 
This is a more general result regarding the ensemble computational nature of 
NMR computing, which allows any spin-direction in the xy-plane to be (in 
theory) resolvable. It is further possible that this alternative embedding could 
be utilised for quantum computations involving only real-valued coefficients. 
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Since such computations are universal [5], it is plausible that this could be 
used to perform ensemble quantum computations in which the state amplitudes 
themselves can be measured du-ectly. 

6 EXPERIMENTAL 

99.8% deuterated CHCI3 was obtained from Aldrich Chemicals. NMR 
spectra of a CHCI3 / H2O mixture were recorded on Bruker Avance II 600 
NMR spectrometer (see Fig. 3), corresponding to a Larmor frequency of 
—600.13 MHz. On-resonant 90° pulse durations were 2.5 ms and Gaussian 
selective pulses at a resolution of 1000 points utilised cut-offs of 1%. All 
spectra were recorded in single scans, allowing recycle delays of at least lis 
between experiments. 
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